Quantum entanglement of a tunneling spin with mechanical 
modes of a torsional resonator 
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Abstract. - We solve Schrodinger equation describing a tunneling macrospin coupled to a tor- 
sional oscillator. Energy spectrum is studied for various quantum regimes. Magnetic susceptibility 
and noise spectrum are computed. We show that entanglement of the spin with mechanical modes 
of a subnanometer oscillator results in the decoherence of spin tunneling. For larger oscillators 
the presence of a tunneling spin can be detected through splitting of the mechanical mode at the 
resonance. Our results apply to experiments with magnetic molecules coupled to nanoresonators. 



Introduction. — There has been enormous progress 
in measurements of individual nanomagnets [T], micro- 
cantilevers and microresonators [2^115). Experiments have 
demonstrated that a mechanical torque induced by the ro- 
tation of the magnetic moment may be used for develop- 
ing high-sensitivity magnetic probes and for actuation of 
micro-electromechanical devices. The underlying physics 
is a direct consequence of the conservation of the total 
angular momentum: spin plus orbital. While this side 
of the effect is clear, the mechanism by which the angu- 
lar momentum of individual spins gets transferred to the 
rotational motion of a body as a whole has been less un- 
derstood. In a macroscopic body it involves complex evo- 
lution of interacting spins and phonons towards thermal 
equilibrium. In that respect the case of a magnetic nano- 
or microresonator is simpler due to the great reduction of 
the number of mechanical degrees of freedom. 

Recently, theoretical study of rotating magnetic 
nanosystems has been conducted within classical [I6HT9] 
and semiclassical [2011^ approaches. When spin is treated 
quantum- mechanically, further reduction of the number of 
degrees of freedom can be achieved in the low energy do- 
main where only the lowest spin doublet that originates 
from the tunneling between spin-up and spin-down states 
is relevant. This would be the case of, e.g., a single- 
molecule magnet. Rigorous quantum-mechanical treat- 
ments have been recently suggested for the problem of a 
tunneling macrospin in a freely rotating body |22j and for 



the problem of a tunneling macrospin coupled to the rota- 
tional modes of a nanoresonator [53] (see Fig. [1]). These 




Fig. 1: System studied in the paper. Macrospin (e.g., a mag- 
netic molecule) is attached to a torsional oscillator such that 
the magnetic anisotropy (quantization) axis is parallel to the 
axis of mechanical rotations. 

two problems have one common feature: The spin tunnel- 
ing becomes suppressed when the body containing the spin 
is too light. The physics behind this effect is quite clear 
Delocalization in the spin space that corresponds to 
tunneling of spin S between spin-up and spin-down states 
reduces the energy by A/2, where A is the splitting of the 
tunneling doublet. Since spin transitions are accompanied 
by the change of the angular momentum they generate 
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rotational motion of the body with the energy / (21) , 
where / is the moment of inertia and L is the mechanical 
angular momentum that is generally of order S. At small 
/ such rotations cost too much energy, so that the tun- 
neling in the ground state should be frozen. This effect 
is conceptually similar to the decoherence and freezing of 
the tunneling of a particle in a double-well potential due 
to dissipation [25] , 

In this paper we solve Schrodinger equation for a 
macrospin coupled to a nanoresonator, Fig. [TJ By consid- 
ering various ranges of parameters of the nanoresonator, 
we reproduce previously obtained results and establish 
connection with the problem of a macrospin in a freely 
rotating body. Qualitatively different behavior is found 
for different ranges of parameters, that can be interpreted 
as a quantum phase transition. We show that the way 
to look for these effects is to study the electromagnetic 
response of the system depicted in Fig. [T] Our most im- 
portant finding is that the coupling of a tunneling spin to 
a mechanical resonator destroys quantum coherence only 
for very small resonators - generally resonators consist- 
ing of just a few atoms. In resonators of greater size the 
coherence is preserved, and the presence of a tunneling 
spin can be detected by observing frequency splitting of 
mechanical oscillations. 

The model. — Consider a model of a tunneling spin 5, 
projected onto the lowest tunneling doublet, in a nanores- 
onator of torsional rigidity k that can rotate around the 
2;-axis [21], see Fig.[TJ 



(heavy particle) that corresponds to the spin tunneling 
between up and down, with the change in the angular mo- 
mentum absorbed by the rotation of the particle. How- 
ever, for a > Qfi the ground state becomes degenerate and 
in the limit a ^ ai (light particle) it approaches J = ztS, 
which means that the spin cannot tunnel. 

In the case of a particle having a restoring force, that 
is the subject of this work, the total angular momentum 
of the spin and the mechanical oscillator is not conserved. 
Conservation of the angular momentum occurs in a larger 
closed system. Still, through the crystal field, tunneling 
of the spin generates mechanical torque acting on the tor- 
sional oscillator [211 [26]. This interaction can seriously 
reduce spin tunneling for both small and large r when the 
oscillator is light. In particular, for small r and large a 
(see below) the ground state is nondegenerate but the gap 
between the ground state and the first excited state be- 
comes very small and the tunneling becomes effectively 
blocked. For this problem Kovalev et al. [23] introduced 
another dimensionless parameter. 
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that is especially useful in the case of large r. One of 
their observations is that at r ^ 1 coupling of the spin 
to quantized rotational vibrations renormalizes the tunnel 
splitting according to 



ioff 



(5) 



2h 



H = 



(1) 

Here Lz = —id^p is the operator of the mechanical angular 
momentum, Iz is the moment of inertia of the resonator, 
= \/k/Iz is the frequency of torsional vibrations, W = 
2SgfiBBz is the energy bias due to the longitudinal field 
Bz, A is the tunnel splitting of spin- up and spin-down 
states due to crystal field, and a are Pauli matrices. As 
we shall see below, the behavior of such a system depends 
on two dimensionless parameters: 
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In the limit of a free particle, r — Q, the total angu- 
lar momentum of the system with respect to the z axis is 
conserved: Jz — Sz + Lz = const. Tunneling of the spin 
changes 5*2 by 25*, and this change is absorbed by the op- 
posite change of Lz- Thus tunneling occurs between two 
quantum states having the same total angular momentum 
eigenvalue J. Computation of the eigenstates of the sys- 
tem reduces to the diagonalization of a 2x2 matrix. The 
resulting spectrum of the system has the ground state with 
J = for [221 



a < ai = 
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(3) 



To solve the quantum mechanical problem of a spin tun- 
neling in a rotating body, it is convenient to use the basis 
that is a direct product of the two-state "up/down" basis 
for the spin and the harmonic oscillator basis for the body. 
Thus we write the system's wave function |^) in the form 



(6) 



m=0 cr=±l 

The coefficients Cma satisfy the Schrodinger equation 

..dCr 

in — - — = 
dt 

n=0 cr' = ±l 

where 

1 
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Bma^mii^aa' 2 Acff 

X (k 

) (8) 



are matrix elements of the Hamiltonian, Eq. ([T]). Here 

E,r,a^nwr{m+\/2)-{l/2)Wa (9) 
are energies in the absence of tunneling and [23) 
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Fig. 2: Distance AE between the ground state and first excited 
state vs Q for different r. 



for m > n and a similar expression with m ^ n for n > m, 



where L„ 



(m— n) 



are a generalized Laguerre polynomials and 

(11) 



A is given by Eq. In particular, 

Kqo = 1, KlO = 1^01 = iA, 



For r ^ 1 only the ground state of the resonator is relevant 
in spin tunneling because energies of all other states are 
too high compared to A In this case one returns to 
a two-state model for the spin with the effective splitting 
([5]). For small r, the spin couples to many oscillator states 
and one has to diagonalize a large matrix. 

Energy spectrum and static susceptibility. 

Setting C„,a ^ C„„Te~'(^/'*)* in Eq. ^ rcsuhs in the 
stationary Schrodinger equation that can be diagonalizcd 
numerically to find energy eigenvalues Ef^. The results for 
the distance /S.E between the ground state and the first 
excited state vs a for different r and W = Q are shown in 
Fig. [2] For r ^ 1 and a > 1 the ground state becomes 
quasidegenerate with very small, although nonzero IS.E. 
This corresponds to the localization of the spin in either 
spin-up or spin-down state. On the other hand, IS.E does 
not exclusively characterize the spin but also contains in- 
formation about the resonator. In particular, for r < 1 
and a — >■ the spin and the resonator effectively decouple 
and Ai? — > tuor, which is the mode of the resonator. On 
the other hand, for r > 1 and a -> one has IS.E — > A, 
which is the spin tunneling mode. 
The spin susceptibility is 



X 



dW 



(12) 



in the limit — )■ 0. For a spin in a massive (non-rotating) 
body one has 



Fig. 3: Reduced inverse susceptibility vs a for different r. 

thus the zero-field susceptibility is xo = V^- For 8- spin in 
a rotating body, the effective splitting, Apfi, can be defined 
through X = 1/Acff, where x = d^a-^jdW follows from 
the exact numerical diagonalization of the Hamiltonian, 



E 

m=0 cr=±l 



(14) 



C^^ma being the coefficients of the wave function corre- 
sponding to the ground state, = 0. The dimensionless 
ratio Xo/x ~ A/A^ff vs a for different r is shown in Fig. 
131 For r ^ 1 and a > 1 the zero-field susceptibility be- 
comes very large because of quasidegeneracy of the "up" 
and "down" spin states. For r ^ 1 Eq. ([5]) is recovered. 



Spin-rotation resonance. — The case vA^TW^ « 
hhJr corresponds to the spin-rotation resonance that leads 
to a strong hybridization of spin and rotational states even 
in the case A ^ 1. In the absence of the interaction be- 
tween the spin and the resonator, A = 0, the lowest four 
energy levels are 



± 



VA2 -I- 



VA2 -f 



(15) 



where the zero-point energy of the resonator has been 
dropped. The hybridized levels are VA2 + 1^2/2 and 
hbjj. — \/A2 -f W"^ jl. The truncated low-energy Hamilto- 
nian matrix has the form 



w_ 

2 



A 







W A-\ A/1 \2 



+ ^ fiA f(l-A2) 



w 

2 







(16) 



V -#*A #(1-A2) 

where Eq. (fTTj) was used. We look for E « 
\/A2 + w hujr/2. Then for A < 1 the equation 
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det(H - M) = simplifies to 
E - i E-hujr 



+ \ A2A2 



UJ± 



E-Eo 



A 

1 ± -- 



A 



2 ^ 



(18) 



This formula provides the splitting of the mechanical and 
spin modes at the resonance. For such a splitting to be 
observable, the quality factor of the mechanical resonator 
must exceed (1+M^^/A^)/A. Eq. (fT8| can also be obtained 
within semiclassical approximation pO] . 

Spin dynamics. — In the problem of a tunneling spin 
embedded in a non-rotating crystal, the parameter A has 
a clear physical meaning of the energy gap between the 
lowest tunneling doublet. When such a spin is prepared in, 
e.g., the spin-up state at t = 0, the probability to find it in 
the same state at an arbitrary moment of time t oscillates 
on t according to (crz(t)crz(O)) = (ctz)* — cos{At/h). When 
the spin is coupled to a light oscillator and r :s> 1, one has 
to replace A Acff. At r ^ 1, coherent probability 
oscillations are destroyed at any non-zero a. 

Spin dynamics is governed by the Schrodinger equation, 
Eq. ([7]), and the time dependence of (ctz) is given by 



(19) 



m=0 a=±l 



(20) 



where Cmu (t) can be expanded over the eigenstates C^; 
as 

the coeficients being determined by the initial condi- 
tion. If at t = the spin was in the "up" state and the 
particle was in its ground state, one has = C*.qi. Com- 
bining these formulas yields the time dependence 
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where 



A., 
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m=0 (T=±l 
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Fourier spectrum of this time dependence, 2|j4^^/|, 
gives the imaginary part of the susceptibility. Plotted vs 
fiw^^i' = E^ — Efj,', it gives an idea of the resonance ab- 
sorption of the rf field by the spin. For r <C 1 the Fourier 
spectrum consists, in general, of many lines. In the limit 
a — ?> the spin and the torsional oscillator decouple; in 
this case only one line of height 1 remains. For a ^ 1 
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(17) 

At the resonance, hur — vA^+W^, the frequencies of the 
transition between the ground state Eq = —y/A^ + W^ /2 0.6 
and the closest excited states become 
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Fig. 4: Fourier spectrum of {(Jz)t 
and 1. 



for r = 0.03 and a = 0, 0.3, 



there is a narrow group of lines with a spread that gives 
rise to spin decoherence due to interaction of the spin with 
the oscillator. At a > 1 decoherence becomes very strong 
and the low-frequency part of the Fourier spectrum corre- 
sponding to mechanical oscillations becomes large. These 
results are shown in Fig. U for r = 0.03. 

At r ^ 1 , there is only one line of height 1 at huj = A^e 
with Acff given by Eq. ([5]) . This is natural because in this 
limit the problem is described by an effective two-state 
model. For r — 1 and small a there is a doublet of lines 
around huj = A because of the resonance interaction 
between the spin and the resonator. 



Discussion. — We have studied energy spectrum, sus- 
ceptibility, and decoherence in a system consisting of 
macrospin rigidly coupled to a torsional mechanical res- 
onator. Our general conclusion is that the coupling does 
not infiuence spin tunneling when the resonator is suffi- 
ciently large and heavy. However, when one approaches 
the atomic size the magneto-mechanical coupling may 
lead to strong decoherence of the spin states. To put 
these statements in perspective, let us consider a magnetic 
molecule of spin 10 embedded in a torsional resonator in 
the shape of a paddle of dimensions 20 x 20 x lOnm'^. As 
in Ref. [23] we shall assume that the paddle is attached 
to the walls by two carbon nanotubes of torsional rigid- 
ity P7] k = 10~"'^^N-m. The moment of inertia of such a 
system is dominated by the paddle, Iz ~ 10~'^^kg-m^, so 

that LUr = y'kjl: 



lO^s"^. The parameter A is then of 
order 10~^. For A/h^ 10^s~^, coherent spin oscillations 
at frequency A/h will not be affected by the coupling to 
the paddle. They will be more likely decohered by the 
coupling to nuclear spins or other environmental degrees 
of freedom in the same manner as for a spin embedded in 
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a macroscopic solid. For A/h > lO^s^^ the parameter r 
will be small. However, a will be very small compared to 
one, and, thus, in accordance with Fig. |4l no decoherence 
of spin oscillations due to coupling with the mechanical os- 
cillations of the paddle will occur in this case either. The 
same will be true even if instead of the paddle one consid- 
ers, e.g., a Mni2 molecule attached to a carbon nanotube 
[28] . The relevant moment of inertia is now that of the 
molecules itself, which for a nanometer size molecule is 
of order 10~^^kg-m^. The corresponding uJr is of order 
lO^^s"^ and A ^ 0.1. The two regimes are now r ^ 1 for 
A/H < lO^^g-i and r < 1, a < IQ-^ for A/H » lO^^g-i 
In both limits the mechanical oscillations should have little 
effect on coherent spin oscillations with frequency A/h. 

The above estimates show that the effects on tunnel 
splitting and spin decoherence arising from spin-rotation 
coupling should not be much of a concern in nanomechan- 
ical setups with large magnetic molecules that have been 
discussed in literature. To have a significant effect on 
the tunnel splitting one should arrive to A > 1. This re- 
quires much smaller moments of inertia, that is, molecules 
much smaller than Mni2. Consider, e.g., a small magnetic 
molecule of spin 10 with the moment of inertia that is one 
hundred times smaller than that of the Mni2 molecule. 
We shall also assume that it is coupled to the walls with 
the torsional rigidity one hundred times smaller than the 
coupling through a carbon nanotube. In this case one still 
gets uir ~ lO^^s"^ but A > 1. Now the regime with r ^ 1 
is achieved at A/?i ~ lO^^s"^, which corresponds to a 1. 
In this case one should expect significant decoherence of 
spin oscillations. The bottom line is that decoherence due 
to a resonator may occur in subnanometer systems but 
it should not be pronounced above the nanometer size. 
This is easy to understand if one notices that Iz^^r in the 
expression A = Sy^2h/IzU!r is the measure of the "macro- 
scopicity" of the resonator. Consequently, IzUt ^ H needed 
to achieve large A generally requires a system of the atomic 
size. For larger resonators, interaction between spin and 
mechanical degrees of freedom reveals itself only near the 
resonance. It results in a very interesting quantum phe- 
nomenon that can be observed in experiment: Splitting of 
the mechanical mode of the resonator containing a tun- 
neling spin. Indeed, in our example with a paddle having 
A ^ 10~^ the splitting of the mechanical mode at the reso- 
nance can be quite significant. For A < Hujr the resonance 
will be achieved at W/h ~ ~ 10^s~^, which for 5* = 10 
corresponds to the magnetic field of order lOG. For, e.g., 
A/H ^ lO^s"^, according to Eq. (fT8| . this will provide 
the splitting in the MHz range that would be possible to 
observe if the quality factor of the resonator exceeds one 
thousand. 
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